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Abstract 

A novel model of intermittency is presented in which the dynamics of the rates of energy transfer 
between successive steps in the energy cascade is described by a hierarchy of stochastic differential 
equations. The probability distribution of velocity increments is calculated explicitly and expressed 
in terms of generalized hyper geometric functions of the type n Fo, which exhibit power-law tails. 
The model predictions are found to be in good agreement with experiments on a low temperature 
gaseous helium jet. It is argued that distributions based on the functions n F$ might be relevant 
also for other physical systems with multiscale dynamics. 
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Our current understanding of fully developed turbulence rests upon two main pillars, 
namely, the energy cascade, whereby energy is transferred from coarser-scaled structures 
to the finer, and the phenomenon of intermittency, which results from the fluctuations of 
the rate of energy transfer. Yet combining these two ingredients into a physically coherent 
model remains an elusive task. In his 1941 theory (K41), Kolmogorov [I] assumed a constant 
rate of energy dissipation, which implies Gaussian statistics for the velocity increments, in 
disagreement with experiments that show heavy-tailed distributions at small scales. The 



lognormal model of intermittency proposed by Kolmogorov [2| in his refined theory based 
on earlier work by Obukhov {3], although in somewhat good agreement with experimental 
data, has been criticized under several grounds . Several other models of intermittency 
have been discussed in the literature 5Ml4|. none of which has been found to be fully 
satisfactory either with respect to their physical basis or in comparison with experimental 
data [15 ]. 



In this paper we present a new model of intermittency where the fluctuating dynamics 
of the rates of energy transfer between successive steps in the energy cascade is described 
by a hierarchy of stochastic differential equations. Under certain reasonable assumptions, 
an integral expression for the stationary probability density function (PDF), p(e r ), of the 
energy flux e r at a given scale r is obtained. From the knowledge of p(e r ), the PDF of the 
velocity increments is then calculated and expressed in closed form in terms of generalized 
hypergeometric functions of the type n Fo, which exhibit power-law tails. The model pre- 
dictions are shown to be in excellent agreement with data from experiments on a turbulent 
gaseous helium jet for several values of the Reynolds number. It is also argued that the 
distributions presented here for the first time are likely to find applications in other physical 
systems with multiscale dynamics. 

According to the energy cascade picture of turbulence, energy is injected at the integral 
scale L and transferred to smaller scales through a hierarchy of eddies of decreasing size, 
until it is dissipated by viscous effects at the Kolmogorov length scale 77. Let us then denote 
by e n the rate of energy (per unit mass) transferred to the scale r = L/b n from the scale 
L/b n ~ l , where b > 1. (Typically, one sets b = 2 but this is not necessary for our analysis.) 
Because the energy flux e n (t) is a fluctuating quantity we seek here to describe its dynamics 
in terms of stochastic processes. On the basis of reasonable physical considerations (see 
below), we propose that the dynamics of e n (t) is governed by the following set of stochastic 
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differential equations (SDE): 



ii = 7j(ej_i - ei) + A^e^t), i = l,...,n, (1) 

where the parameters 7« and ki are assumed to be constant in time and £i(£) are mutually 
independent white noises. The quantity €q appearing in Eq. ([1]) for i — 1 represents the rate 
of energy fed into the system at the integral scale L and is considered fixed. 

The terms in the right-hand side of Eq. (CQ) have a clear physical interpretation. For 
instance, the deterministic term represents the unidirectional coupling between successive 
steps of the energy cascade. Owing to this coupling, if we were to neglect the fluctuating 
term in Eq. (TjQ) then all quantities e, would relax to the constant value eo, thus recovering 
the K41 theory. By the same token, Eq. ([1]) implies that the average energy flux is scale 
independent, in the sense that in the stationary regime one has (e,) = eo for all i. The choice 
of the noise term is also a natural one since we expect a multiplicative noise in a cascade 
process. This ensures, in particular, that if the quantities 6i(t) are initially positive then they 
remain nonnegative for all times. To see this, note that if 6i(t) were ever to become negative 
it would have to cross zero, since it is a continuous process. But if = at some time, 
then Eq. (TjQ) implies that q > and so e» will be 'reflected' back to the positive range. (Of 
course, the rate of energy transfer e r cannot assume negative values if it is to be identified 
with the local average rate of energy dissipation, as first suggested by Obukhov 

The model defined in Eq. ([I]) bears some resemblance to shell models of energy cascade 
in turbulence 16|, where one seeks to describe the energy- cascade mechanism by a set of 
coupled nonlinear ordinary differential equations that are consistent with the Navier-Stokes 
(NS) equation. Our model is more of a phenomenological nature in that it incorporates 
the fluctuations of the rates of energy dissipation explicitly via a set of coupled stochastic 
differential equations. We note, however, that it is possible to give a heuristic derivation 
of the deterministic term in Eq. jT]) from the scale- by-scale energy budget equation [3| , if one 
assumes localness of the energy transfer. In the same vein, the noise term can be justified 
from symmetry considerations and from the positivity requirement on (see above). A 
related approach based on energy-balance equations was used in to obtain a Langevin 
description of the energy of eddies of different sizes, but here the resulting SDE's are highly 
nonlinear. Our model, in comparison, is written in terms of the energy transfer rate, is 
linear, and has the further advantage that it yields an analytical expression for the PDF of 
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the velocity increments which is in excellent agreement with experimental data, as we will 
see shortly. 

Understood in the sense of the Ito stochastic calculus, Eq. flTJ constitutes a set of n linear 
SDEs that can be solved exactly 20(. Such an approach, however, is not very useful for us 
here since it is not easy to obtain the stationary joint probability distribution p{e\, ...,e n ) 
from this exact solution. Considering the stationary Fokker-Planck equation for p(e\, e n ) 
is not very helpful either, since this equation cannot be easily solved. Thus, an alternative 
approach is needed to compute the stationary PDF for e n . Here we will take advantage of 
;he separation of the characteristic time scales at the different steps of the energy cascade 



181 ] . To be specific, we make the following assumption: 7" 1 <C ln-i "C • • • <C <C 7] -1 - 
On the basis of this hypothesis, we can derive the stationary PDF for e n from our dynamical 
model, as follows. 

Consider first Eq. (]T|) for i = n. Since the dynamics of e n has a characteristic time much 
shorter than that of e n _i, it is reasonable to assume that before e n _i has time to change 
appreciably the flux e n relaxes to a quasi-stationary regime described by a conditional PDF, 
p(e n |e n _i), obtained assuming e n _i fixed. In other words, the marginal distribution for e n 
can be written as a superposition of distributions p(e n |e„_i) with different values of e n -\. 
p(e n ) = / °°p(e n |e n _i)p(e n -_i)c/e n _i. Implementing this procedure recursively up to the first 
step of the energy cascade, we obtain 

/*oo /*oo n 

p(e n ) = / ■••/ Y\p{e i \e i - 1 )de 1 - ■ -dtn-x. (2) 
Jo Jo i=1 

The distribution p(ej|ej_i) can be obtained by solving the stationary Fokker-Planck equation 
associated with Eq. (jrj, holding ej_i fixed. This yields an inverse-gamma distribution 

P(*l*-i)= m + 1) ^ e s , (3) 

where 

A = 2 7l /kl (4) 

With the knowledge of the PDF of the energy flux e n , we can now derive the PDF of the 
longitudinal velocity increments, 5 r u = u(x + r) — u(x), at a given scale r. To this end, we 
express the marginal distribution for 5 r u as 

P(5 r u) = J p(e r )P(6 r u\e r )de r , (5) 



where P(5 r u\e r ) is the conditional probability distribution of 8 r u for a fixed value of e r . 
Since intermittency stems from the fluctuations of e r , it is reasonable to assume that the 
statistics of the velocity increments for fixed e r is described by a Gaussian distribution. This 



assumption is supported by experiments 

P(S r u\e r ) 
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221 ] . We then write 
(5 r u) 2 ' 



V2 



exp 



7T(J- 



2a 2 



(6) 



where a 2 is the (random) variance of 5 r u for a fixed va' 



a 2 with the local average energy dissipation rate e r [9j. More formally, however, we write 



ue of e r . In general, one can associate 



a 2 = ((6 r uf\e r ) = ((6 r u) 2 ) 



(7) 



where the second identity is to be understood in the measure-theoretic sense, meaning that 
the random variable (Su 2 \e r ) is a coarser version of 5u 2 23 1. 

From Eqs. (J5j) — it then follows that the PDF of 5 r u normalized to unit variance can 
be written as 



P(6 r u) 



= exp 



(5 r u 



r.\2 



2e r 



de r 



(8) 



where 5 r u = 5 r u/ \J ((5 r u) 2 } is the normalized velocity increment and e r = e r /eo is the 
normalized energy flux. Because we assume a Gaussian of zero mean in Eq. ([6]), our model 
describes only the symmetrical part of the PDF of the velocity increments, whose non- 
Gaussianity is a signature of intermittency 24|. The asymmetry (skewness) of the PDFs 



is thought to be connected with vortex folding and stretching [24j and is (for the moment) 
left out of the model. The idea expressed in Eq. (JSJ) of writing the PDF of the small-scale 
velocity fluctuations as a mixture of large-scale (Gaussian) distributions has been used by 
several authors with different weighting distributions, such as the gamma distribution [9|, the 



lognormal distribution 25H27j. and the chi-square distribution 28]. A related description 



based on a Fokker-Planck equation for the conditional distribution of velocity increments 
was introduced in 



29(. In comparision to these previous works, the novelty of our approach 
is that we model the dynamics of the energy fluxes and then derive (rather than postulate) 
its distribution, from which the PDF of velocity increments can be obtained explicitly, as 
shown next. 

Upon substituting Eqs. (T5]) and (j3J) into Eq. (JS}, and performing a sequence of changes of 



variables, one can show 



171 ] that the resulting multidimensional integral can be expressed 
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in terms of known higher transcendental functions: 

r(A + 3/2) 



P(5 r u) 



n 



nF {ax,...,a n ;- SrU ), (9) 

ZP! • • • p n 



.1=1 

where a, = & + 3/2 and n F (ai, a n ; — x) is the generalyzed hypergeometric function of 



order (n, 0) 30]. The first two members of the family n Fo yield elementary functions, namely, 
o-F is the exponential function and iF is related to the so-called g-exponential: iF (l/(l — 
q),x) = exp q (x), where exp (x) = [1 + (1 — q)x} 1 ^ 1 ~ q \ We thus see that the distributions 
in Eq. (Q give a rather natural generalization of the Gaussian (n = 0) and the g-Gaussian 
(n = 1) distributions. One important property of the function n Fo, for n > 0, is that it has 
an asymptotic expansion of the form n F (ai, a„; — x) oc Yu7=i c i x ~ ai (1 + 0(l/x)), 
as x —>■ oo. Thus, the distributions P(x) above comprise a general class of power-law tail 
distributions with finite variance. (This seems to be the first time that distributions based 
on the functions n Fo, with n > 1, appear in the literature.) 

Returning to our intermittency model given in Eq. (JTJ, we now make the simplifying 
assumption that the parameters /?, are the same throughout the cascade: fti = (3. This 
implies, in particular, that the distribution p(ej|ej_i) given in Eq. ([3]) is scale invariant in 
the sense that it has the same functional form regardless of the cascade level. In view of the 
discussion in the preceding paragraph, it then follows that P(8u) in this case has a single 
power-law tail: P(5u) ~ Su~^ 2 ^ +3 \ for Su 3> 1. 

Next we compare the model with experimental velocity measurements on the axis of 
a low temperatu re g aseous helium jet. For details about the experiments the reader is 



referred to Refs. 
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331 ] . From the recorded data sets, typically with 10 7 points each, we 
computed the velocity differences 5u between two consecutive measurements. In Fig. [TJ we 
show the (symmetrized) histograms of velocity increments for four values of the Taylor-scale 
Reynolds number, namely, R\ = 463, 703, 885, and 929, together with the corresponding 
PDFs (solid lines) predicted by our model. The agreement between the theoretical curves 
and the experimental data in Fig. [1] is remarkable. 

In Fig. [T] we chose n as the smallest value necessary to fit satisfactorily the data, in the 
sense that increasing n gives no further improvement of the fit. This is illustrated in Fig. [2] 
where we show the experimental histogram for R\ = 703, together with the theoretical fits 
for n = 1,2, 3, 4, 5. Here for each n we chose the value of (3 so as to fit the largest possible 
range of the data. As n increases, the agreement between the theoretical curve and the 
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FIG. 1: (color online). Histograms of velocity increments (circles) computed from measurements 
on the axis of a gaseous helium jet for four different Reynolds numbers: R\ = 463, 703, 885, 
929 (from bottom to top). The corresponding solid lines are the theoretical PDFs for n = 4 and 
f3 = 6.6, 6.5, 6.4, 6.2. The curves have been arbitrarily shifted in the vertical direction for clarity. 

data improves considerably, up to a point where further increasing n makes no practical 
difference. Note, in particular, that the g-Gaussian (n = 1) is in a rather poor agreement 
with the experimental data, failing most notably to fit the tails 34]. The data shown in 
Fig. [2] is for the smallest separation resolved by the experiment, namely, r = 6.5 /im, as 



obtained from Taylor's 



case is rj = 3.4 fim 



rozen turbulence hypothesis. Since the Kolmogorov scale in this 
321 ]. this indicates that our model is apparently valid down to the 
intermediate dissipation range 18] . We have verified that the model is also able to fit the 
PDFs of velocity increments computed at larger separations, with the number of cascade 
steps required to fit the data decreasing when r increases, as expected. (More details will 
be given elsewhere \y\-) 

As a final point, let us briefly consider the structure functions predicted by our model. 
Using the known properties of the inverse-gamma distribution, one can show that 



(e p n ) 



p-1 

n 







(10) 



Recalling that r = L/b n , it then follows that the moments of e r naturally obey a scaling 
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du / [<du 2 >] m 
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FIG. 2: (color online). Distribution of velocity increments for R\ = 703 (circles) and corresponding 
fits (solid lines) for different values of the number n of steps considered in the energy cascade. The 
curves (from top to bottom) are for n = 1 and /3 = 1.9 (red), n = 2 and (3 = 3.4 (blue), n = 3 and 
f3 = 5.0 (green), n = 4 and /3 = 6.6 (black), and n = 5 and (3 = 7.9 (magenta). 



relation 

<^>oc(£f P , (11) 

where t p = — log 6 (/3//3 — i). Note, however, that the velocity structure functions do 
not necessarily exhibit scaling, since {(S r u) 2p ) = (2p — l)\\((5 r u) 2 ) p {e p ) and we impose no a 
priori scaling for the second moment of the velocity increments. (A more detailed discussion 
about the scaling properties of our model will be left for a forthcoming publication 
We also note in passing that our intermittency model recovers the lognormal model in the 
limit of an infinite cascade. Indeed, if we take the limits n — > oo and (3 — > oo, in such way 
that of = n/(3 remains finite, then Eq. ( TTUj) becomes 



(e p r ) 



e e ) 



(12) 



which are precisely the moments of a lognormal distribution lnA/"(eo, of) 35j. Our model 
thus provides a dynamical context where the lognormal model naturally arises. 

In conclusion, we have presented a new cascade model of intermittency in fully developed 
turbulence based on a hierarchy of stochastic differential equation for the energy fluxes at 
different scales in the cascade. The model is derived from a physically reasonable set of 
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assumptions and produces an analytical formula for the PDF of the velocity increments in 
terms of the generalized hypergeometric functions n F , which fits extremely well the exper- 
imental data. We conjecture that distributions based on n F are likely to find applications 
in other systems whose dynamics entails multiple spatial or temporal scales, such as frag- 
mentation processes, biosystems, and financial data. 
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